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The intrinsic shape of galaxy bulges 


J. Mendez-Abreu 


Abstract The knowledge of the intrinsic three-dimensional (3D) structure of galaxy 
components provides crucial information about the physical processes driving their 
formation and evolution. In this paper I discuss the main developments and results in 
the quest to better understand the 3D shape of galaxy bulges. I start by establishing 
the basic geometrical description of the problem. Our understanding of the intrinsic 
shape of elliptical galaxies and galaxy discs is then presented in a historical context, 
in order to place the role that the 3D structure of bulges play in the broader picture 
of galaxy evolution. Our current view on the 3D shape of the Milky Way bulge and 
future prospects in the field are also depicted. 


1 Introduction and overview 

Galaxies are three-dimensional (3D) structures moving under the dictates of gravity 
in a 3D Universe. From our position on the Earth, astronomers have only the op¬ 
portunity to observe their properties projected onto a two-dimensional (2D) plane, 
usually called the plane of the sky. Since we can neither circumnavigate galaxies nor 
wait until they spin around, our knowledge of the intrinsic shape of galaxies is still 
limited, relying on sensible, but sometimes not accurate, physical and geometrical 
hypotheses. 

Despite the obvious difficulties inherent to measure the intrinsic 3D shape of 
galaxies, it is doubtless that it keeps an invaluable piece of information about their 
formation and evolution. In fact, astronomers have acknowledged this since galaxies 
were established to be island universes and the topic has produced an outstanding 
amount of literature during the last century. 
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In this paper I discuss the main developments and results in the quest to bet¬ 
ter understand the 3D shape of galaxy bulges. Given the limited space avail¬ 
able in this chapter, I have not elaborated on the concept and definition of a 
bulge, leaving this discussion to another chapter in this volume. In the same 
way, I have deliberately not included the intrinsic shape of boxy/peanut (B/P) 
structures located in the centre of disc galaxies wh ich some authors associate 
to galaxy bulges ( Liitticke. Dettmar & Pohletl 200r)h . Currently it is well estab¬ 
lished that these structur es are actually part of the bar and intimately re lated 
to their secular evolution ( Combes & Sanded 198ll: Chung & Bureau . 2004 ). As 
bars evolve, stars can be moved perpendicular to the disc plane due to a coher¬ 


ent bending of the bar producing its characte ristic shape (IDebattista et al.L 12004 


Martinez-Valpuesta, Shlosman & Helleil 200^. B/P structures share the same ph* 


tomet ric and kinematic properties of bars (IMendez-Abreu et al.LuOOSbclErwin & Debattista , 


120131). 


On the other hand, I have included a historical review of the evolution of our 
knowledge of the intrinsic shape of elliptical galaxies. The properties of elliptical 
galaxies and thos e of intermediate/massive galax y bulges have been often consid¬ 
ered to be similar ( Wvse. Gilmore & Franxi 1997 ). This is particularly true when re¬ 
ferring to their surface-brightness distributions and shapes. Indeed, it has been com¬ 
mon in the literature to rely on both simulations and observati ons of elliptical galax¬ 
ies to interpret the observational properties of bulges ('e.g.. iKormendv & Bender . 


2012 ). 


This paper is structured as follows. In Sectional describe the basic geometric 
considerations of the problem and set up the notation used throughout the chapter. 
In Section[3]l review our current knowledge on the intrinsic shape of both elliptical 
and disc galaxies. Section |4] introduces the advantages and drawbacks of studying 
galaxy bulges with respect to ellipticals and a historical perspective of their 3D 
shape measurements. In Section|5]l summarize the evolution of the concept of the 
Milky Way bulge and its intrinsic 3D shape. Section |6] addresses the importance 
of numerical simulations to understand the physical processes that shape galaxy 
ellipsoids. Finally, in Section|2]I sketch out the current view on the intrinsic shape 
of bulges and explore future prospects. 


2 Setting up the Scene 

This section briefly summarizes the basic notation and geometrical considerations 
to be used during this chapter. 

Let ix,y,z) be the Cartesian coordinates on the reference system of the galaxy 
with the origin in the galaxy centre, the x—axis and y—axis corresponding to the 
principal equatorial axes of the ellipsoidal component, and the z—axis correspond¬ 
ing to the polar axis. Therefore, if A, B, and C are the intrinsic lengths of the ellipsoid 
semi-axes, the corresponding equation of the bulge on its own reference system is 
given by 
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Fig. 1 Schematic three-dimensional view of the ellipsoid geometry. The bulge ellipsoid, the disc 
plane, and the sky plane are shown in red, blue, and green, respectively. The reference systems 
of both the ellipsoid and the observer as well as the LON are plotted with thin solid lines, thin 
dashed lines, and a thick solid line, respectively. The bulge ellipsoid is shown as seen from an 
arbitrary viewing angle (left panel), alo ng the LOS (central panel), and along the polar axis (i.e., 
the z—axis; right panel). Extracted from llMendez-Abreu et alll201(J) . Reproduced with permission 
from Astronomy & Astrophysics, © ESO. 
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Let (x',y,z') now be the Cartesian coordinates on the observer reference system. It 
has its origin in the galaxy centre, the polar z'—axis is along the line of sight (LOS) 
and points toward the galaxy, (x',/) represents the plane of the sky. 

The equatorial plane (x,y) of the ellipsoid and the plane of the sky [x! ,y') inter¬ 
sect in the so-called line of nodes (LON). The angle between both planes, i.e., the 
angle subtended between z and ^ is defined as the inclination 0 of the ellipsoid. The 
remaining two Euler angles which allow for the transformation from the reference 
system of the galaxy to that of the sky are defined as: i) (j) is the angle subtended 
between the x—axis and the LON in the ellipsoid equatorial plane, and ii) y/ is the 
angle subtended between the y—axis and the LON in the plane of the sky. It is often 
useful to choose the x'—axis to be along the LON, consequently it holds that y/ = 0 
(see Figure[T]i. 

It is well known t hat the projection of a triaxial ellipsoid onto the plane of the sky _ 

descr ibes an ellipse ( Contopoulos[ll956HStarkLll977l : lBinnev[ll985l : rFranx. Illingworth & de Zeeuw . 


19911) . which is usually written as 
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+ yl 

2 + ^2 


= 1 


( 2 ) 


where Xe and represent the axes of symmetry of the projected ellipse, a and 
b are the corresponding semi-major and semi-minor axes of the ellipse. The ob¬ 
served ellipticity of the ellipse can be easily derived from the apparent axis ra¬ 
tio as e = 1 — bja. The Xe axis forms an angle 5 with the LON (twist angle). 
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which for convenience is usually made to correspond with the x' -axis. It is worth 
noting that both the apparent axis ratio {q = bjd) and the orientation of the el¬ 
lipses (5) depend only, and unambiguously, on the direction of the LOS, i.e., on 
9, 0, and i/r, and on the intrinsic shape of the ellipsoid, i.e., A,B, and C, see 
Simonneau, Varela & Munoz-Tunon ( 1998h for the full derivation. 


Based on this simple geometric representation, if we assume a galaxy is com¬ 
posed of a set of triaxial emitting ellipsoidal shells, which are concentric and coax¬ 
ial (same axes of symmetry) but non-homologous (intrinsic semi-axes vary with the 
distance to the centre), their projections onto the plane of the sky are concentric el¬ 
lipses, but non-homologous and non-coaxial. Therefore, the twisting of the galaxy 
isophotes can be explained just as an effect of the projection of non-homologous 
triaxial ellipsoids ( Williams & Schwarzschild , 19791) . 


3 Historical background on the intrinsic shape of galaxies 


Elliptical galaxies are structurally the simplest stellar systems where mathematical 
techniques can be applied to recover their intrinsic 3D shape. Thus, the huge amount 
of literature on the subject is not surprising. In fact, the continuously increasing 
availability of better measurements of the apparent axis ratios of elliptical galax¬ 
ies have motivated great debate over the years. On the other hand, the similarities 
betw een the photorn etric properties of intermediate/massive bulges and ellipticals 
(e.g., lOadottil l2009l) have usually motivated an extrapolation of the results on the 
intrinsic 3D shape of ellipticals and their implications on galaxy formation and evo¬ 
lution onto the bulges of disc galaxies. In this section I revisit our current knowledge 
on the intrinsic shape of elliptical galaxies (Sect. 13.1b and, for the sake of complete¬ 
ness, of disc galaxies (Sect. 13.2b to put in context the historical background on the 
intrinsic shape of bulges. 


3.1 Intrinsic shape of elliptical galaxies 


3.1.1 Photometric approach 


The first attem pt to derive the intrinsic shape of elliptical galaxies was done by 
Hubble ( 1926h . At that time, it was already realized the importance of relying on 
statistical methods to recover the 3D shape of galaxies. In fact, Hubble obtained the 
frequency of intrinsic short-to-long axis ratio under the assumption that elliptical 
galaxies were oblate ellipsoids with random orientations with respect to the LOS. 

Since then, this statistical approach based on the measurement of the apparent 
axis ratio distribution (AARD) and the assumption that the 3D intrinsic shape is an 
ellipsoids of revolution, either oblate or prolate, has been extensively used in the 
literature. For the sake of clarity I briefly outline here the basic statistical concepts. 
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Let us assume the basic geometry proposed in Sect.l^and define both the intrinsic 
ellipticity, Q = BjA, and intrinsic flattening, F = C/A, of the ellipsoid as the corre¬ 
sponding intrinsic axis ratios in the (x,y) and {x,z) planes, respectively. Therefore, 
in the case of either a pure oblate (Q = 1) or pure prolate (Q — F) ellipsoid the Eq. 
[T]can be described by one single parameter. If the polar axis of the ellipsoid forms 
an angle (9) with respect to the LOS then the apparent axis ratio of the projected 
ellipse can be written as 


F^ sin^ 0 -f cos^ 0 = 


if oblate 
if prolate 


(3) 


Under the realistic assumption of randomly distributed orientations and using Eq. 
[3]where q = q{9), the probability P{q\F)dq that a galaxy with intrinsic axis ratio F 
is observed with an apparent axis ratio in the range {q^qF dq) is 


P{q\F)dq 


sin 0 dq 
\dq/d9\ 


(4) 


At this point, the AARD {q), can be related to the intrinsic probability distribu¬ 
tion ^ (E) by 


C P{q\P)^{F)dF (5) 

JO 

The relation between the known (observed) frequency of galaxies of apparent 
axis ratio l^{q) to the unknown frequency ^ [F] of galaxies with intrinsic axis ratio 
F can be written such as 


C(?) = 


■L 


‘I ^ {F)dF 


L 


V'(1-E2)(^2_^2) 

? ^{F)F^dF 


V'(1-E2)(^2_^2) 


if oblate 


if prolate 


(6) 


Bas ed on this approach and using the hypothesis of oblateness. ISandage, Ereeman & Stokes 
( 1970l) derived the intrinsic distribution of flattening ^(E) for different Hubble 
types ranging from ellipticals to Sc. They found that the observed axis ratios of 
168 elliptical galaxies present in the Ref erence Catalog of Bright Galaxies (RCl) 

( de Vaucouleurs & de Vaucouleurs , 19641) were well reproduced using a skewed bi¬ 
nomial distribution of oblate ellipsoids given by 




E-Eo 


exp[-a{F -Fo)] 


(7) 


wi th main parame ters Eq = 0.58 and j3 = 0.31 (Eigure|2] left panels). 

BinnevI ( 1978 ) used the same sample but introducing the prolate approach. 


Adopting the same functional form for ^ (E) he found values of Eq = 0.40 and 
j3 = 0.71. However, even if using arbitrary analytical representations of ^{F) 
can turn out in a good fit of the AARD, in principle they do not have a physi- 
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cal motivation. This approximation was improved by Noerdling^ ( 1979f) by solv - 
ing Eq. |6]using the non-parametric inversion technique proposed by ILucvI (119741) . 
His results show how under the hy pothesis of oblateness the ^{F) distribution of 
ISandage, Freeman & Stokes ( 19701) was correct, but he also noticed that a prolate 
distribution peaking at around F ^ 0.7 would produce a good representation of the 
data as well. 

At the same time, some kinematic findings led to the suggestion that the struc¬ 
ture of elliptical galaxies can be represented by neither oblate nor prolate ellipsoids 
of revolution. In fact, the low ratio between rotational velocity and veloc i ty disper¬ 
sion f ound in flat systems ( Bertola & CanaccioliL 1975 ; Illingworth , 1977t PetersonL 
1978h or the rotation mea sured along the minor axis of some elliptical galaxies 


( Schechter & Gunni 19791) were interpreted as resulting from a triaxial structure. 
From the photometric point of view, the t wisting of the inner isophotes of elliptical 
galaxies was kno wn since the early work oflEvansI ( l^M ) and it was lat ter confirmed 
in several works ( Fillet]^ 1960^ C^tei^ 1^78 ; Bertola & Gallettai 1^79^ _ 

As a consequence, Benacchio & Galletta ( 1980l) and Binnev & de Vaucouleurs 
(1198 ih showed that the AARD could be satisfactorily accounted for also in terms of 
a distribution of triaxial ellipsoids. Nevertheless, these works still presented signifi¬ 
cant differences in the predicted number of spherical galaxies mainly due to the dif¬ 
ferences in the original samples. Other groups reached similar conclusions a nalysing 
higher quality data coming from new CCD detectors ( Fasano & Viol 1991 ). 

A new step forwa rd in the methodolo gy to recover the intrinsic 3D shape of 
galaxies was done by Fall & FrenkI ( 1983h . They showed how the inversion of the 
integral equations for oblate and prolate ellipsoids (Eq.|6l) can be performed analyt¬ 
ically, resulting in 


^(F) = -v/1-F2. 


1 cF 

F Jo 


‘idq ^ 

dq 


1 rF qdq djq^Q 


fo 




dq 


if oblate 


if prolate 


(8) 


Using this analytical i nversion and the largest samp l e of g alaxies to that date 
(2135 elliptical galaxies), Fambas, Maddox & Foved^ ( 1992 ) demonstrated how 
neither oblate nor prolate models could adequately reproduce the data. Contrarily, 
triaxial ellipsoids with intrinsic axis ratios selected from ID Gaussians provided an 
adequate fit to the dat a. They found a best fit with Q = 0.95 and F = 0.55. A similar 
approach was used bv iRvdenI (Il992h on a smaller sample of 171 elliptical galaxies. 
She used a 2D Gaussian combining both intrinsic axi s ratios obtaining Q = 0.98 
and F — 0.69. The same sample was later analysed by Tremblay & Merritt ( 1995h 
using a non-parametric technique to test the triaxial hypothesis. They confirmed 
previous results that discarded a distribution of intrinsic shapes compatible with 
axisymmetric ell ipsoids thus f avouring triaxial distributions. Similar conclusions 
were reached by Rvden ( 19961) on a larger sample using the same non-parametric 
approach. 

During these years it became increasingly clear that the distribut ion of intrinsic 
flattenings of elliptical galaxies was broad and possibly bimodal ( Fasano & 
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199 it Rvdenl 1992t Tremblay & MOTi^ 1995 : Rvden . 1996 ). In fact, combining 
the galaxy sample des cribed in R^den|j 19921) w it h a new sample of brightest cluster 
galaxies (BCGs) from Lauer & PostmanI ( 1994 ). Tremblay & Merritt ( 19961) found 
that the AARD of galaxies brighter than Mg — —20 was different from that of the 
less luminous ones. This reflected a difference in the shape of low-luminosity and 
high-luminosity ellipticals: fainter ellipticals are moderat ely flattened and ob late, 
while brighter ellipticals are rounder and triaxial. Recently, Fasano et al. ( 2010l) also 
found that eyen if both normal ellipticals and BCGs are triaxial, the latter tend to 
haye a more prolate shape, and the tendency to prolateness is mainly driyen by the 
central dominant (cD) galaxies present in their sample. 

The next qualitatiye leap in studies of the intrinsic shape of elliptical galaxies 
happened with the adyent of the Sloan Digital Sky Suryey (SDSS). With respect 
to preyious statistical analyses, SDSS improyed not only the number of galaxies 
under study (an order of magnitude larger) but also the quality and homogeneity 
of the photometry. All these improyements allowed to study the dependence of the 
intrinsic shape with other galaxy properties such as the lu minosity, colour, physi¬ 
cal size, and enyironmen t. Using data from the SDSS-DR3 (lAbazaiian et al.L 120051) 


Vincent & RydenI ( 2005 ) found that bright galaxies (Mr < —21.84) with a de Vau- 


couleurs profile haye an AARD consistent with a triaxiality parameter in the range 
0.4 <T < 0.8, where T = (1 — Q^)/{\ —F^), and mean flattening 0.66 <F < 0.69. 
The faintest de Vaucouleurs galaxies are best fit with prolate ellipsoids (T = 1^ 
with mean flattening F = 0.51. Using the SDSS-DR5 (lAdelman-McCarthy et al 


2007ft . iKimm & Yil (l2007h were able to reproduce the AARD by using a combina¬ 


tion of oblate, prolate, and triaxial galaxy populations. Following the early work of 
Tremblay & Merritt ( 19961) . they assumed each population haying a Gaussian dis¬ 


tribution of their intrinsic axis ratios. The best fit to the AARD was found using a 
fraction of O:P:T=0.29:0.26:0.45 (Oblate:Prolate:Triaxial ) with a best triax i al dis- 
tribution with axis rat ios Q = 0.92 and F = 0.78. In 20 08, Padilla & StraussI ( 2008h 
used the SDSS-DR6 ( Adelman-McCarthy et all 2008h to deriye the intrinsic shape 
of ellipticals with the main improyement of taking into account the effects of dust 
extinction. They found that the AARD of elliptical galaxies shows no dependence 
on colour, suggesting that dust extinction is not important for this sample. The full 
population of elliptical galaxies was well characterized by a Gaussian distribution 
in the equatorial ellipticity with mean Q = 0.89 and a lognormal distribution of 
the flattening wit h mean F = 0.43, which corresponds to slig htly oblate ellipsoids 
in agre ement with Vincent & RydenI ( 20051). In a rece nt paper, Rodriguez & Padilla 
( 20131) haye used the SDS S-DR8 jA ih^a et all 201 ih and the morphological infor¬ 
mation from Galaxy Zoo dUintott et al. . 2011 ) finding that elliptical galaxies haye a 
mean yalue of F = 0.58 (Figured right panels). They concluded that the increase 
in F is mainly due to the remoyal of the spiral galaxy contamination thanks to the 
Galaxy Zoo morphologies. A historical summary in tabular form of all these mea¬ 
surements is shown in Table [T] 

Owing to the ill-posed problem of deriying the 3D intrinsic shape of elliptical 
galaxies, its historical perspectiye is mainly weighted toward statistical methods. 
As preyiously showed in this section, the inyentiyeness of astronomers, the deyel- 












































































J. Mendez-Abreu 


opment of statistical methods, and the advent of large surveys have significantly 
improved our knowledge of the intrinsic shape of elliptical galaxies. Other meth¬ 
ods based on the photometric study of individual galaxies have also been devel¬ 
oped but to a smaller extent. One of the pioneering works to derive the intrinsic 
shape of an individual ellipti cal using its observed ellipticity and isophotal twist 
was done by Williams ( 198lh . They modelled the elliptical galaxy NGC 0523 as¬ 
suming a given intrinsic density distribution and finding that the preferred mod¬ 
els were prolate in the external regions but increasingly mixed (oblate and prolate) 
towards the centre. This idea was further d eveloped by other authors using more 
comp lex model s of the density distribution (Fasano . 1995 : Thakur & Chakrabortvl 
200 ih . In 2008, Chakrabortv. Singh & Gaff^2008 ) estimated the shapes of 10 el¬ 


liptical galaxies with apparent ellipticities e < 0.3, finding that radial differences 
in the triaxiality parameter can be tightl y constrained to values 0.29 < AT < 0.54. 


Chakrabortv. Diwakar & PandevI ( 201ll) extended this analysis to 3 very flat galax¬ 


ies with ellipticity e ^ 0.3 or more. They found values of the intrinsic flattening of 
these galaxies around F ^ 0.5. 


3.1.2 Kinematic approach 


Determining the distribution of the 3D intrinsic shape of elliptical galaxies is also 
possible by co mbining photometric and kinematic information. In a first attempt, 
Binnevl ( 'l985h used simple kinematical models to understand the ratio of rotational 
motion along both the major and minor isophotal axes of the galaxy. Using a sample 
of 10 ellipticals he found that ellipti cal galaxies were not well repre s ented by ax- 
isymmetric oblate or prolate models. Franx. Illingworth & de Zeeuw ( 199 ih revis¬ 
ited this approach by using a larger sample of 38 elliptical galaxies and studying the 
probability distribution of photometric ellipticities and kinematics misalignments. 
In particular, they explored the possibility that the angular momentum could not be 
aligned with the polar axis of the galaxy but it may have any orientation within the 
plane containing the short and the long axis (x,z). They found that a variety of mod¬ 
els was able to reproduce the observations. Models with all galaxies being triaxial 
with well-aligned angular momentum were indistinguishable from models with all 
galaxies being oblate with nonaligned angular momentum. 

A different standpoint to statistical studies implies an investigation into the in¬ 
trinsic shape of elliptical galaxies using det ailed individual dynamical modelling 
of the galaxy kinematics. iTenies et akl (119931) modelled the photometric and stellar 
kinematic measurements of three elliptical galaxies adopting a specific form for the 
intrinsic density and streaming motions. They found tightly constrained geometries 
with 0.7 < Q < 0.8 and 0.4 < F < 0-6. This methodology was further improved in a 
series of papers by Statler ( Statleil 1994at Statler & Frvl 19941: State, 1994b ). He 
showed how using not only their apparent shapes and velocity field misalignments, 
but also the velocity field asymmetry, it is possible t o place tighter constr aints on 
the intrinsic shape of ellipticals. Using this approach Bak & Statler ( 2000l) derived 
the intrinsic shape of 13 elliptical galaxies finding that although photometric studies 
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Fig. 2 Composite figure showing the evolution of the galaxy samples used in the derivation of 
the intrinsic shape of ellipticals and discs. Upper panels: histograms of the AARD for ellipticals 
and spiral galaxies. The overplotted curves are predicted ratios for various assumptions of the 
distribution of intrinsic flattening. On the right, the assumed intrinsic d istribu tion corresponding 
to the curves on the left. Extracted from ISandage. Freeman & Stok^ il970h . Reproduced with 
permission, © AAS. Bottom panels: best fit models to AARD compared t o the observations. Top : 
spirals. Bottom: ellipticals. Left: galaxies selected only by fracDeV, see lAbazaiian et al.l l l2005h 
for definition. Right: galax ies selected by Galaxy Zoo morphology and fracDeV. Extracted from 
iRodrfguez & Padilla ^20131) . Reproduced by permission of Oxford University Press. 




give similar results for the flattening, none is able to put real constraints on triaxiality 
even when large samples are studied, hence demonstrating the need to include kine¬ 
matic data in the models. Figure|3]show the probability distribution of intrinsic axis 
ratio for 9 galaxies with significant rotation in their sample. It is clear that most of 
the galaxies can be well described by nearl y oblate models but some of th e m pre sent 
significant triaxiality or even prolateness. Ivan den Bosch & van de VenI (l2009l) in¬ 
vestigated how well the intrinsic shape of elliptical galaxies can be recovered by 
fitting realistic triaxial dynamical models to simulated photometric and kinematic 
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Table 1 Historical summary of the intrinsic shapes of elliptical galaxies. 


Year 

(1) 

N. Galaxies 
(2) 

Hypothesis 

(3) 

Q 

(4) 

F 

(5) 

Reference 

(6) 

1970 

168 

Oblate 

1 

0.58 

[1] 

1978 

168 

Prolate 

0.4 

0.4 

[2] 

1979 

168 

Oblate/Prolate 

1/0.7 

0.55/0.7 

[3] 

1980 

348 

Triaxial 

0.81 

0.62 

[4] 

1981 

196 

Oblate/Prolate/Tri axial 

1/0.62/0.79 

0.62/0.62/0.57 

[5] 

1992 

2135 

Triaxial 

0.95 

0.55 

[6] 

1992 

171 

Triaxial 

0.98 

0.69 

[7] 

2005 

26994 

Triaxial 

0.66-0.85 

0.66-0.69 

[8] 

2007 

3922 

Oblate/Prolate/Tri axial 

1/0.72/0.92 

0.44/0.72/0.78 

[9] 

2008 

303390 

Triaxial 

0.89 

0.38 

[10] 

2013 

112100 

Triaxial 

0.88 

0.58 

[11] 


Notes. (1) Year of publication of the paper. (2) Number of elliptical galaxies in each sam¬ 
ple. (3) Hypothesis used to derive the intrinsic shape of the ellipticals. (4) Mean value 


of the intrinsic ellipticity. (5) Mean value of the intrinsic flattening. (6) Reference of 
the corresDonding naner: 111 ISandage. Freeman & Stoke.sl 1197(1), 121 iBinnevI Il978t). 131 

iNoerdlinged lll979l). [41 iBenacchio & Gallettal ||198( 

). [51 iBinnev & de VaucouleursI il981h. 

161 Lambas. Maddox & Lovedavl 1 19921). 17] RvdenI 

19921). 181 Vincent & RvdenI l2005l). 191 

iKimm&Yil 120071). llOllPadilla & StraussI l2008l). 1111 

Rodriguez & Padillal J2013h. 


observations. They found that for axisymmetric galaxies, the models are able to 
exclude triaxiality but the intrinsic flattening is nearly unconstrained. On the other 
hand, the shape of triaxial galaxies can be accurately determined when additional 
photometric and kinematic complexity, such as the presence of isophotal twist or a 
kinematicall y decoupled core is obs erved. 

Recently, IWeiimans et al.l (120141) studie d the intrinsic shape o f the early-type 


galaxies described in the ATLAS3D survey (ICappellari et al.Ll2011h . Using a purely 


photometric approach and assuming axisymmetry, they found that the fast rotator 
population was much flatter than the slow rotator population, as expected from their 
dynamical status. Moreover, when the kinematic misalignment is included as a con¬ 
straint in the analysis, they demonstrated that fast rotators are still better represented 
to oblate ellipsoids. 


3.2 Intrinsic shape of disc galaxies 

In this section I briefly summarise our current understanding about the intrinsic 3D 
shape of discs. Bulges are embedded into the disc light and axisymmetry is usu¬ 
ally a requirement to derive the bulge intrinsic shape. However, although the discs 
of lenticular and spiral galaxies are often considered to be infinitesimally thin and 
perfectly circular, their intrinsic shape is better approximated by flattened triaxial 
ellipsoids. 
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NGC 4261 
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Fig. 3 Posterior probability densities in the plane of intrinsic triaxiality, T, and flattening, cl (F i n 
this chapter), for each of the nine galaxies that show significant rotation in iBak & StatleJ i2000h . 
Contours indicate the 68% and 95% highest posterior density regions. In each panel, round prolate 
galaxies are at the top left, flattened oblate galaxies at bottom right, and objects in between are tri- 
axial. Most galaxies are well represented by oblate models but prolate and tr iaxial are also allowe d 
in many galaxies, e.g., NGC 741, NGC 4486, or NGC 7626. Extracted from iBak & Statlej l l2000h . 
Reproduced with permission, (c) AAS. 


The disc flattening, deflned analogously as for ellipticals (Sect. 13.1b . can be 
directly determined from edge-on galaxies. It depends both on the wavelength at 
which discs are observed and on gala xy morphological type. Indeed, galactic discs 


becom e thicker at longer wavelengths jPalcanton & BernsteinLl2002 ; Mitronovaet^ 
2004) and late - type s pirals have thicker discs than early-type spirals ( Bottinelli et all 


1983HGuthrielll992[) . 


Determining the distribution of both the intrinsic flattening and ellipticity of 
discs is possible by a statistical analysis of t he AARD of randomly oriented spi- 
ral galaxies. Similarly for elliptical galaxies, Sandage. Freeman & Stokei ( 1970t) 
analysed the spiral galaxies listed in the RCl. They concluded that discs are cir¬ 
cular with a mean flattening of {F) = 0.25. However, the lack of nearly circular 
spiral galaxies (o' ~ 1) rules out that discs have a perfectly axisymmetric shape. In- 
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deed. 


Binggeli ( 1980l) . Benacchio & Gallettal ( 1980h . and Binnev & de VaucouleursI 


(Il981 ) have shown that discs are slightly elliptical with a mean intrinsic ellipticity 


{\ — Q) =0.1. These early findings were based on the analysis of photographic 
plates of a few hundreds of galaxies. They were later confirmed by measuring 
ellipticities of several thousands of o bjects in CCD images and digital sc ans of 
plates obtained in wide-held surveys. iLambas. Maddox & Lovedavl (119921) found 
that pure oblate models failed to reproduce the AARD of spiral galaxies, whereas 
nearly oblate mode ls with F ^ 0.2 and Q ^ 0. 9 produ ce a good ht with values 
similar to those of Sandage, Freeman & Stokes ( 1970l). These values were con- 


hrmed later on by different authors ( Fasano et al. . 1993 : Alam & Rvdenl 2002; 


Rvdenl 2004). Like the battening, the intrinsic ellipticity depends on the morpho¬ 


logical type and wavelength. The discs of early-type spirals are more elliptical than 
those o f late-type spir als and their median ellipticity increases with observed wave¬ 
length ( Rvdenl 20061) . Furthermore, luminous spiral galaxies tend to have thic ker 
and rounder discs than low-luminosity spira l galax ies IPadilla & StraussI (l2008h . In 


Sanchez-Janssen. Mendez-Abreu & Aguerri ( 2010l) they studied the role of stellar 


mass in shaping the thickness of galaxy discs. They found that the intrinsic thick¬ 
ness distribution of discs has a characteristic U-shape and identify a limiting mass 
M* ft! 2 X lO^MfT) below which low-m ass galaxies start to be systematically thicker. 
Recently, Rodriguez & Padilla ( 20131) analyse a sample of 92923 spiral galaxies ex¬ 
tracted from the SDSS-DR8, and taking into account the effects of dust in their 
analysis, they found a distribution of battening with mean F = 0.27 and ellipticity 
Q = 0.22, i.e., disc are less round than in previous studies (Figured right panels). 

Despite the large effort made to understand the intrinsic 3D shape of galaxy discs, 
it is still unclear whether the inferred slight triaxiality could be due to the presence 
of substructure in galaxy discs or if it really rebects truly triaxial potential in spirals. 


4 The intrinsic shape of extragalactic bulges. 


The study of the intrinsic shape of bulges presents similarities, advantages, and 
drawbacks with respect to that of elliptical galaxies. Bulges are ellipsoidal sys¬ 
tems located in the centre of disc galaxies, thus, the main drawback with respect 
to elliptical galaxies is that their analysis requires the isolation of their light distri¬ 
butions from other structural galaxy components. However, it is worth noting that 
a similar problem is faced in elliptical galaxies w hen debning a charac teristic ra¬ 
dius to measure the global axis ratio of the galaxy (IFasano & Violll99ll) . The most 
common approach to identifying a global axis ratio for the bulge is by perform¬ 
ing a photometric decomposition of the galaxy surface-brightness distribution. In 
this method, the galaxy light is usually modelled as the sum of the contributions 
from the different structural c omponents, i .e., bu l ge and disc, and eventually lenses . 


bars, spiral arms, and rings (IPrieto et al.L 1200It iLaurikainen, Salo & ButaL l2005h . 


A number of two-dimensional parame tric decompos ition techni ques have been de - 
veloped to this aim, such as; GIM2D ( Simard . 1998 ). GALFIT ( Peng et al. . 2002), 
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BUDD A ( de Souza . Gadotti_&dos_Aiijos|^ 200^, G ASP2D ( Mendez-Abreuetal 


2008ah . GALPHAT dYoon, Weinberg & Katzl boTH) . or IMFIT (lErwinLl2015h . On 


the other hand, the main drawback on the study of galaxy bulges, i.e., the presence 
of other components such as the main disc, represents in turn the main advantage. 
The presence of the galactic disc allows for accurately constraining the inclination 
of the galaxy. Hence, under the assumption that the two components share the same 
polar axis (i.e., the equatorial plane of the disc coincides with that of the bulge) it 
allows for the determination of the inclination of the bulge. This is crucial to solve 
one of the main concerns when dealing with elliptical galaxies. 


4.1 Photometric approach 


Galaxy bulges were initially thought as axisymmetric ellipsoids placed at the cen¬ 
tre of dis c galaxies . The first piece of photometric evidence against this idea was 
given by LindbladI (1956). He showed a misalignment between the major axes of 
the disc and bulge in M31, realising that this would be impossible if both the disc 
and bulge were oblate. This photometric misalignment is similar to the isophote 
twist obs erved in elliptical galaxies and us ed as an indication of triaxiality in these 
system s ( Williams & Schwarzschildl 1979h . The extensive study undergone bv iKenj 
( 19841) showed that the twisting isophotes between the central and outer parts of disc 
galaxies are quite common, but it was not until 1986 when Iz^itskv & O (Il986l) 
properly studied the deviations from axisymmetry in the bulges of spiral galaxies. 
They found bulge-to-disc misalignments in their sample of 11 spiral galaxies hence 
confirming the high incidence of non-axisymmet ric bulges in ordinary spirals and 
placing some parallelisms with elliptical galaxies. Beckman et al. ( 199 ih also found 
compelling photometric eyidence for triaxiality in the bulge of NGC 4736. 

The first quantitatiye estimation of the intrinsic 3D shape of galaxy bulge s using 
a statistical approach was performed by Bertola. Vietri & Zeilingei ( 199 ih . They 
measured the bulge AARD and the misalignments between the major axes of the 
bulge and disc in a sample 32 SO-Sb galaxies. Under the hypothesis that discs are 
circular, they found that these bulges are triaxial with mean axial ratios {Q) = 0.86 
and {F) = 0.65. Interestingly, they also demonstrated that a random projection of 
the probability distribution function of the bulge s axis ratios fit sufficiently well t o 
the AARD of the elliptical galaxies presented in iBinney & de VaucouleursI (Il98ll) . 
T he results were in t erpret ed as both populations of objects haying the same origin. 
Fathi & Peletierl ( 20031) deriyed the intrinsic ellipticity of bulges by analysing the 


deprojected apparent axis ratio of the galaxy isophotes within the bulge radius. This 
work did not assume any geometrical model for the galaxy but only that the disc be 
circular. They found {Q) = 0.79 and {Q) = 0.71 for the bulges of 35 early-type and 
35 late-type disc galaxies, respectiyely. Despite the dif ferent methodologies, these 
results were in good agreement with preyious results by lBertola, Vietri & Zeilinger 
( 199 ih . Along the same lines, non e of the 21 disc galaxies with morphol ogical types 
between SO and Sab studied by iNoordermeer & yan der HulstI (l2007h harbours a 
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truly spherical bulge. They reach this conclusion by assuming bulges to be oblate 
ellipsoids and comparing the isophotal axis ratio in the bulge-dominated region to 
that measured in the disc-dominated region. A mean fla ttening (F) = 0.55 was ob¬ 
tained which is slightly lower than the value found by iBertola. Vietri & Zeilingei 
(Il99lh . 

The number o f galaxy bulges under study i ncreased by an order of magnitude 
with the work of iMendez- Abreu et alJ ( 2008ah . They measured the structural pa¬ 
rameters of bulges and discs of a sample of 148 early-to-intermediate spiral galaxies 
using a 2D photometric decomposition. They computed the probability distribution 
function of the intrinsic ellipticity from the bulges AARD, disc ellipticities, and 
misalignments between bulges and discs position angles. They suggested that about 
80% of the sample bulges are triaxial ellipsoids with a mean axial ratio (B/A) = 
0.85, confirming that bulges are slightly triaxial structures. 

The vertical extension of gala xy bulges remains usually hidden from obse rva- 
tions except for edge-on galaxies. Mosenkov, Sotnikova & Reshetnikov ( 2010h ob¬ 
tained a median value of the flattening (F) = 0.63 for a sample of both early- and 
late-type edge-on galaxie s using near infrared photometry. T hese results match well 
with the early findings bv iBertola. Vietri & Zeilingen(ll99lh . 

As well as for elliptical galaxies a number of works have attempted to quantify 
the intri nsic shape of individual bulges using only pho tometric data. The pioneering 
work of Varela. Munoz-Tunon & Simmoneau ( 19961) used a combination of geo¬ 
metrical deprojection and photometric inversion to work out the actual shape of the 
galaxy bulge in NGC 2841. They found that a family of triaxial ellipsoids with vari¬ 
able axis ratios is necessary to explain the photo metric properties of its bulge. In 
1998, Simonneau. Varela & Munoz-TunorJ ( 19981) derived a set of equations defin¬ 
ing the three intrinsic axes of a triaxial ellipsoid as a function of the measured ge¬ 
ometry of a galaxy bulge and disc (axis ratios and position angles) and the unknown 
Euler angle 0 (see Sect.lHfo r definition). This seminal paper promoted the work of 


Mendez-Abreu et al.l (l201(]h . They introduced a new method to derive the intrinsic 


shape of bulges based upon the analytical relations between the observed and in- 
trinsic shapes of bulges and their surroundi ng discs. Using the equations derived in 
Simonneau. Varela & Munoz-Tunon ( 1998h and introducing physical constraints on 
the accessible viewing angles, they found the following relation between the intrin¬ 
sic semi-axes of the bulge and their observed properties 


2 sin (2(^c) p2 ^ ^2^^ - (j)g)^ (1 _|_ g 2 

re 


(9) 


where 0 b, 0c and Fg are functions of the observed qu antities a, b, 5, and 0, see 
equations 12, 13, and 43 of Mendez-Abreu et al. ( 2010l) . Therefore, Eq. |9] directly 
relates the intrinsic 3D shape of the bulge with its observed properties. Unfortu¬ 
nately, the relation between the intrinsic and projected variables also depends on 
the spatial position of the bulge with respect to the disc on its own reference sys¬ 
tem (i.e., on the 0 angle) and therefore, as well as for ellipticals, a deterministic 
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so lution of the problem canno t be given. However, the statistical analysis provided 
in Mendez-Abreu et al. ( 2010f) allows us to obtain the probability distribution func¬ 


tion of both semi-axis ratios, Q and F, for every single bulge, thus imposing tight 
constrain ts on its actual shape. Applyin g this technique to the sample of bulges pre¬ 
sented in iMendez- Abreu et alJ (l2008al) they found a bimodal distribution of the tri- 
axiality parameter (FigurelH left panel). In particular, bulges with Sersic index n<2 
exhibit a larger fraction of oblate axisymmetric (or nearly axisymmetric) bulges, a 
smaller fraction of triaxial bulges, and fewer prolate axisymmetric (or nearly ax¬ 
isymmetric) bulges with respect to bulges with n>2. Despite no correlations being 
found between the intrinsic shape of bulges and other properties such as bulge lu¬ 
minosity or velocity dispersion, the differences with the bulge surface-brightness 
di stribution hint towards the pres ence of different bulge populations as suggested 
by Kormendv & Kennicuttl ( 200^ . 




Fig. 4 Composite figure showing the similar himodal distrihution of triaxiality parameters from 
observations (left panel) and simulations (right panel). Left panel: di stribution of the triaxiality 
parameter T obtained from the sample of iMendez- Abreu et all 1201 Oh (continuous line) and for 
a simulated sample with both 30% an d 100% of bulges hosting a nuclear bar (dashed and dot¬ 
ted lines), respectively. Extracted from iMendez-Abreu et all l2010h . Reproduced with permission 
from Astronomy & Astrophysics, (c) ESO. Right panel: distribution of both dissipation al (hatched 
histog ram) and dissipationless (solid line) mergers remnant triaxiality parameter from ICox et all 
I 2 OO 6 II . In both panels oblate galaxi es have T = 0, p rolate galaxies have T = I, and all values in 
between are triaxial. Extracted from ICox et all l200dl . Reproduced with permission, (g) AAS. 
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4.2 Evidences of triaxiality from kinematic measurements 


Early kinematic st udies of galaxy bulges were sho wn to rotate more rapidly than 


elliptical galaxies ( Kormendv & IllingworthL [l982 1. In fact, the kinematic proper- 


which are flattened by rotation with little or no anisotropy ( 

Davies & Illingworth. 

1983 

; Jarvis & Freeman, 1985; Fillmorel 19861; Corsini et al.. 

19991; Pignatelli et al.. 

2001 

in a 

). However, there are also kinematic evidences supporting a triaxial shape 
non-negligible fraction of these bulges. In 1989, two independent works of 


same conclusion about the triaxial bulge of the Sa galaxy NGC 4845. Using a com¬ 
bination of photometric and kinematic measurements they restrict the intrinsic axis 
ratio of its bulge to 2 = 0.74 and F = 0.6. Their works were mainly supported by the 
presence of non-circular gas-motions in the galaxy centre. In a non-axisymmetric 
potential, the shape of the rotation curve will depend on the position of the LOS 
and the major axis of the non-axisymmetric component. A slowly rising rotation 
curve or one in whic h a bump of extreme velocities is seen near the centre are indi¬ 
cations of triaxiality ( Gerhard. Vietri & Kend 19 89h. Based on the se considerations. 


and building on the early work of Lindbladp 1956h . BermanI ( 2001 ) demonstrated the 


presence of a triaxial bulge in the Andromeda galaxy (M31) by using a hydrodynam- 
ical simulation to match the observed properties of the gala xy. Further evidences 


for no n-cir cular gas motion in g alaxy centres can be found in iFalcon-Barroso et alJ 


( 2006h and Pizzella et al. (2008). Other kinematic evidence for the existence of tri¬ 


axial buIges_comes from the presence of velocity gradients along the galaxy minor 


axis. 


Corsini et al. ( 2003h found minor axis rot a tion in 80% of their early-type spi¬ 


ral sample. In a series of papers, Coccato et al. ( 2004 2005 1 found that 60% of the 
unbarred galaxies show a remarkable gas velocity gradient along their optical minor 
axis. This was achieved by combining their own data with t hat present in the liter¬ 
ature (Revised Shapley-Ames Catalog of Bright Galaxies) (ISandage & Tammann , 
198lh . 

Despite the importance of adding kinematic information to determine th e intrin¬ 
sic sha pe of the bulges, and contrary to the works on elliptical galaxies (e.g.- IStatler . 
1994al) . there is not a well-established methodology to quantify the degree of triax¬ 
iality of bulges using the combined photometric and kinematic information, yet. 


4.3 Polar bulges 


Polar bulges, as well as their analogous polar rings ( Whitmore et al. , 1990h . are elon¬ 
gated structures perpendicular to the plane of the galaxy disc. A common signature 
of both the orthogonally decoupled bulge systems and the polar ring galaxies is that 
both contain a structural component whose angular momentum vector is roughly 
parallel to the major axis of the host galaxy. 
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Vertical elongation is not a common feature of bulges. Indeed, most bulges 
can be assumed to be flattened by rotation (see Sect. 14.2b . Furthermore, orthog¬ 
onally decoupled bulges are usually not even allowe d in most statistical works 
since th e condition A > B > C is co mmonly used, see iBertola, Vietri & Zeilinger 
( 199 ih . Mendez-Abreu et 51 ( 2010h relaxed this condition and found that only 
18% of the observed bulges have a probability > 50% of being elongated along 
the polar a xis with no bulges reaching a p robability > 90% . In fact, to date 
NGC 4698 (lBertolaet allll999h . NGC 4672 dSarzi et all lioooll . and UGC 10043 
( Matthevys & de Griisl 20041) are the only spiral galaxies known to host a prominent 
bulge sticking out from the plane of the disc. 

The case of NGC 4698 is particularly intriguing since it hosts also a polar nuclear 
stellar disc aligned with its polar bul ge and thus p e rpend icular to the main disc. 
This galaxy was recently revisited by Corsinietal. (201 ^ and its intrinsic shape 
was derived using the methodology proposed by Mendez-Abreu et al. ( 2010h . They 
found a slightly triaxial polar bulge elongated along the vertical direction with axis 
ratios Q = 0. 95 and F = 1 . 60. Th is result agrees well with the observed kinematics 
presented in Bertola et alJ ( 1999ll and with a model where the nuclear disc is the 
end result of the acquisition of external gas by the pre-existing triaxial bulge on the 
principal plane perpendicular to its shortest axis and perpendicular to the main disc 
of the galaxy. 


5 The Intrinsic Shape of the Milky Way Bulge 


Owing to its vicinity, the Galactic bulge has always been targeted as the ideal bench¬ 
mark for structure, kinematic, and stellar populations studies of bulges. In fact, it can 
be studied at a unique level of detail, in comparison to external galaxies, thanks to 
the possibility of measuring the properties of individual stars. However, our inside 
view of the Galaxy generally restricts our knowledge to pencil beam areas around the 
Galactic centre due to either the high extinction, the crowding, or the superposition 
of multiple structures along the LOS, making studies of the inner Galactic regions 
challenging. The structure of the Galaxy has accounted for a significant amount of 
literature in the past and the topic has come back in the limelight in recent years. In 
this section I briefly review the Galactic bulge topic focusing on its intrinsic shape 
heading the readers to other chapters in this volume for more information about its 
stellar content and kinematics. 

_ In recent decades it has become clear that the Galaxy is a barred system ( Blitz & SpergelL 

1991 ; Lopez-Corredoira, Cabrera-Lavers & Gerhard! 2005 ) and that most likely its 
central regions are dom inated by a boxy bulge created by vertical instabilities 


within the Galac tic bar ( Dwek et al. , 1995 ; Martinez-Valpuesta & Gerhardl 2011 


Ness et alil2013h . The historical evolution of our knowledge of the intrinsic struc¬ 


ture of the Galactic bulge has been written by a succession of progressively larger 
scale, deeper sensitivity photometric and spectroscopic surveys. 
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The first attempt to underst and the shape of the Galactic bulge was made by 


de Vaucouleurs & Pencel (119781) . They found that models ranging from spherical 


to F = 0.6 were able to represent well both the distribution of globular clus- 
ters around the Gala ctic centre and the infrared isophotes observed at 2.4/rm 


dMaihara et all 119781) . The flattening of the Galactic bulge was then further con- 


strai ned with the arrival of the In frare d Astronomical Satellite (IRAS) . Using IRAS 
data, Harmon & Gilmor j ( 1988 ) and Whitelock. Feast & Catchpole (1991) found 
values of the int rinsic flattening spanning 0.6 < T < 0.8 using JHK near-infrared 
bands. Similarly, Kent, Dame & Fa^ ( 199lh found that, at first order, the Galactic 
bulge can be represented by an oblate ellipsoid with F = 0.61 using data from the 
Infrared Telescope (IRT). _ 


Th e picture changed drastically with the advent of the COBE satellite (IHauser et al 
1990l) . The new striking image of the Mi lky Way (Figure!^ prov ided by the D IRBE 

elRli 


experi ment on board of COBE allowed iBlitz & Snergell (1199 Ih . and later on [Blitz 
(Il993h . to find the first direct evidence for a bar at the Galactic centre. Interest¬ 
ingly, they also found the presence of a triaxial bulge structurally distinct from the 
main bar. The modelling of this triaxial bulge was performed by different teams 
with different sets of data in the subsequent years. C onsequently, differe nt axis ra- 
tios represented as 1:Q:F w ere found: 1:0.33:0. 22 (lDweket al.[ 1995 ). 1:0.6:0.4 
(iBinnev. Gerhard & Snergell Il997h 1:0.43:0.29 (iStanek et al.l Il997h . 1:0.38:0.26 
(Freudenreich, 1998) , 1:0.5 4:0.33 ( Lonez-C orredoiraet al., 2 000l). l:('0.3-0.4):0.3 
( Bissantz & GerhardLl2002h . 1:0.5:0.4 (iLopez-Corredoira. Cabrera-Lavers & Gerhard . 
200.5h . In general, these values implied the Galactic bulge to be a triaxial structure 


with a tendency to prolateness, thus not in agreement with the triaxial/oblate picture 
outlined in Section|4|for extragalactic bulges. 

Although the i dea of a triaxial bulge worked w ell at f i rst ord er, the boxy shape 
noticed earl i er by Kent. Dame & Fa^ ( 199lh and Kent! ( 1992 ) and confirmed by 
Dwek et ^ ( 1995 ) was not recovered by a triaxial ellipsoid. In the meanwhile, dif¬ 


ferent scenarios came up to explain these differences and accou nt for t he con tinu- 
ously increasing kinematic and stellar populations information. lAlardI (1200 ll) sug¬ 
gested the presence of two different bars in the Galaxy by analy sing data from the 
Two Micron All Sky Surve y (2MASS) ( Skrutsk ie et al.L 20061) . Another possible 
scenario was worked out bv iBabusiaux et al.l (1201 Ol) suggesting a model composed 
by a classical bulge i n the centre and a boxy bulge in the outer parts. 

Shen et al. ( 2010l) proposed a simple model yet backed up by the high qual- 


ity stellar kinem atics provided by the Bulge Radial Velocity Assay (BRAVA) 
( Rich et al. . 2007h . Using N-body simulations they found no evidence for a classical 
bulge in the Galaxy but the bulge appears to be only part of the bar and therefore not 
a separated component. Figure|6]shows that the i nclusion of a class ical bulge greatly 
worsens the model fit to the data. Models from Shen et al. ( 201(]h rule out that the 
Milky Way has a signifi cant classical bulge with mass >15% o f the disc mass. 

(201T) demonstrated how 
(l2007h agrees with a sce- 


Following this line, iMartinez -Valpues ta & Gerhard 
the star counts measurements by ICabrera-Lavers et al 


nario composed by a single bar and a boxy bulge. More recent measurements of 


star counts from the VISTA Variables in The Via Lactea (VVV) ((Gonzalez et al 
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Fig. 5 False-colour image of the near-infrared sky as seen by the DIRBE. Data at 1.25. 2.2, and 
3.5 /tm wavelengths are represented respectively as blue, green and red colours. The image is 
presented in Galactic coordinates, with the plane of the Milky Way Galaxy horizontal across the 
middle and the Galactic centre at the centre. Credits: E. L. Wright (UCLA), The COBE Project, 
DIRBE, NASA. 


201 ih . 


metallicity gradient s from the Abun dances and Radial velocity Galactic 
Origins Survey (ARGOS) ( Ness et al. . 2013), or stellar kinematics from BRAVA 


have also been reconciled within this p icture dGerhard & Martinez-Valpuestall2012t 


have also been reconciled witnin this p i 

iMartihez-Valpuesta & Gerhardl 20131) . 


6 The 3D shape of bulges in numerical simulations. 


The ihtrihsic shape of bulges keeps importaut ihformatioh about their formatioh his¬ 
tory, with differeht merger, accretioh aud assembly sceuarios resultihg ih differeht 
shapes. Hehce, the comparisoh of measured ihtrihsic shapes with the output from 
humerical simulatiohs represeuts au ihtrihsic way to gaih ihsights oh their forma¬ 
tioh. However, humerical resolutioh problems have ofteu hampered these studies 
aud our ihterpretatioh of the shapes of bulges is usually restricted to the aualysis of 
si mulated ell i ptical galaxies. 

Cox et alJ (120061) studied the structure of ellipsoidal remuauts formed by either 


major (equal-mass) dissipatiohless or dissipatioual mergers of disc galaxies. They 
fouhd a bimodal distributioh of the triaxiality parameter ih their remhaht ellipti¬ 
cals (see right pahel ih Figure |4|. Thus, dissipatiohless remhahts are triaxial with 
a tehdehcy to be more prolate ahd with a meah triaxiality parameter T = 0.55, 
whereas dissipatiohal remhahts are triaxial ahd tehd to be much closer to oblate 
with triaxiality T = 0.28. This simulated bimodal distributioh was compared by 
Mehdez-Abreu et al. ( 2010l) to the triaxiality measured ih their sample of 115 galaxy 
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Fig. 6 Best models fits to the BRAVA stellar kinematics using different hypothesis on the classi¬ 
cal bulge mass. Mean velocity (top panels) a nd velocity dispers ion (lower panels) profiles of all 
available kinematic observations presented in IShen et al.l ll2010h . The left two panels are for the 
Galactic latitude b = 4 strip; the middle two panels are for the b = 8 ; and the right two panels 
are for the 1 = 0 minor axis. The heavy black lines represent the model without a classical bulge. 
The red, green, and blue lines are for models whose classical bulges have masses of 8%, 15%, and 
30%, respectively, of the disk mass. Including a classical bul ge significantly wo rsens the model fits 
to the data, especially along the minor axis. Extracted from Shen et al] j2010h . Reproduced with 
permission, AAS. 


bulges (FigurelHi. They concluded that both major dissipational and dissipationless 
mergers are requir ed to explain the v ariety of shapes found for bulges. The detailed 
study presented bv ICox et al. ( 20061) is c onsistent with previous studies of dissipa¬ 
tionle ss and dissipational me rgers (e.g., Barnes, 1992^^ ^rnauistL 1992t Springe], 
2000l) . However, the study of Gonzalez-Garcfa & Balcells (2005) they found how 


the degree of triaxiality of the elliptical remnants in dissipationless mergers also de¬ 
pends on the morphology of the progenitor spirals. The presence of central bulges 
on the progenitor galaxies produce remnants which tend to be more oblate whereas 
bulgeless progenitors lead to highly triaxial remnants which seems inconsistent with 
observations. Therefore, the comparison between simulations and observations are 
still subject to the range of initial conditions explored by numerical simulations. 

On the other hand, even if the similarities between bulges and ellipticals have 
prompted observers to compare the measured properties of bulges to the prop¬ 
erties of simulated elliptical galaxies, the formation path of bulges is likely a 
more co mplex process involving the in t eraction with other g alaxy structural com¬ 
ponent^ ^^mendx^^^nnicutl 2004 : Athanassoulai 200.51) . The recent work by 


Tapia et al.l (120141) has started to fill the gap on studies about the intrinsic shape of 


galaxy bulges from numerical simulations. They analysed a set of A— body Simula- 
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Fig. 7 Intrinsic shape of bulges and elliptical galaxies obtained from numerical simulations. A 
comparison with observed bulges is shown in the second panel. The blue and green stars in all 
panels represent the bulge remnants after suffering intermediate/minor mergers. The location of the 
progenitor bulges is shown with orange stars. The elliptical remnants of major mergers with pure 
exponential stellar discs (black circles) and containing 40% of gas (red circles) are also shown. First 
panel: intrinsic ellipticity b (Q in this chapter) versus the intrinsic flattening c (F in this chapter) 
Secon d panel: as panel 1 but adding the observed distribution of bulges in iMendez- Abreu et alj 
^201011 (black diamonds). Third and fourth panels: triaxiality p arameter as a function of the intrinsic 
ellipticity and flattening. Extracted from iTapia et d] j2014ll . Reproduced with permission from 
Astronomy & Astrophysics, (gi ESO. 


lions of intermediate and minor dry mergers onto SOs to understand the structural 
and kinematic evolution induced by the encounters. In their experiments, the pro¬ 
genitor bulges are nearly spherical. The remnant bulges remain spherical as well 
{Qr^F > 0.9), but exhibiting a wide range of triaxialities (0.20 < T < 1.00), re¬ 
marking how the definition of this shape parameter is too sensitive to nearly spher¬ 
ical systems. Figure |7] (second panel) shows how the axis ratios derived from these 
simulations (open stars) are hard l y reco ncilable with the observations (black di¬ 
amonds) by iMendez-Abreu et^ (l2010l) . Still, the strong triaxiality agrees with 
the structure of e lliptical remnants resulting from major-to-intermediate mergers 
( Cox et all 2006h . 


7 Concluding remarks and future prospects 

I present here a review of our current understanding of the intrinsic 3D shape of 
galaxy bulges. The approach taken in this review is largely observational and fol¬ 
lows the historical development of the field. Thus, a journey through the past and 
present of our knowledge on the intrinsic shape of other galaxy ellipsoids such as 
elliptical galaxies or galaxy discs was needed to put the problem in context. The 
major conclusions of this review are: 

• The observational data representing the whole population of elliptical galaxies is 
consistent with a mixed model, combining partly oblate and partly prolate galax¬ 
ies, although a more likely alternative point towards at least some fraction of the 
ellipticals being triaxial ellipsoids. Triaxiality is also supported by several photo¬ 
metric and kinematics properties, as well as for detailed modelling of individual 
galaxies. 
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• The intrinsic shape of ellipticals shows a dependence on galaxy luminosity. 
Bright ellipticals are in general triaxial with a tendency to be rounder whereas 
faint ellipticals are more flattened with a tendency to be oblate ellipsoids. 

• Even if uncertainties due to the lack of number statistics have been overcome 
with the advent of recent surveys, the data can still be reproduced by a wide 
variety of intrinsic shape distributions. Furthermore, a proper interpretation of the 
data is complicated by the fact that the AARD and kinematic misalignments are 
often a function of the radius. Therefore it is generally impossible to characterize 
the full shape of a single elliptical galaxy with only one or two parameters. 

• Galaxy discs are, in general, well represented by nearly oblate models with Q ~ 
0.9. Their intrinsic flattening is also well constrained to values spanning 0.2 < 
F <03. 

• The population of galaxy bulges can be modelled as slightly triaxial ellipsoids 
with a tendency to be oblate. This population has typical intrinsic flattenings of 
F ~ 0.65. However, individual galaxies can have a variety of intrinsic flattenings 
with some extreme cases sticking out the plane of the disc, these are called polar 
bulges. 

• The distribution of the triaxiality parameter of galaxy bulges is strongly bimodal. 
This bimodality is driven by bulges with Sersic index ii > 2. According to nu¬ 
merical simulations they can be explained assuming a combination of major dis- 
sipational and dissipationless mergers during their formation. 

• Despite previous findings showing a triaxial bulge in the Milky Way, more recent 
studies have found that is more likely a boxy bulge produced by the vertical insta¬ 
bilities of the Galactic bar. Owing to recent kinematic measurements a classical 
bulge with mass > 15% of the disc mass can be ruled out. 

Despite the study of the intrinsic shape of elliptical galaxies has a long track 
record, our knowledge of the 3D shape of bulges is still in its infancy. Therefore, 
further work on the topic is needed to fully exploit its possibilities. A few guidelines 
to this future prospects are outlined in the following: 


From a photometric point of view, even if new methodologies have been devel¬ 
oped they need to be applied to larger samples of galaxy bulges. The number of 
elliptical galaxies recently analysed to recover their intrinsic shape is several or¬ 
ders of magnitude larger than the current samples of galaxy bulges. Farge number 
statistics have led to the discovery of important relations for ellipticals galaxies, 
such as the different shapes of bright and faint ellipticals, and similar studies 
can be crucial for galaxy bulges. This is particularly relevant in the current pic¬ 
ture of bulge formation with a different population of classical and pseudobulges 
dependent of the galaxy mass ( Fisher & Drorvl 2011 ). 

An even more promising path, already explored in elliptical galaxies, is the use 
of combined information from photometric and kinematic data. In particular, the 
common use of integral field spectroscopy is now providing an exquisite detail 
of the stellar and gaseous kinematics on large sample of galaxies. This wealth of 
information together with the development of galaxy dynamical modelling can 
provide a proper understanding of the intrinsic shape of galaxy bulges. 
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It is doubtless that the comparison of the derived intrinsic shape of bulges with 
the state-of-the-art numerical simulations is a promising way to gain insights on 
the formation and evolution of bulges. However, there is still a lack of simulations 
with a large variety of initial and physical conditions interested on a structural 
analysis of the different galaxy components, and in particular, in the intrinsic 
shape evolution of galaxy bulges. 

Historically, galaxy bulges were thought as single-component objects at the cen¬ 
tre of galaxies. This picture is now questioned since different bulge types with 
different formation paths ha ve been found coexisting within the same galaxy (see 
Mendez-Abreu et al. . 20l4 and references therein). A proper separation of dif¬ 


ferent bulges types, as well as the identification of possible unresolved nuclear 
structures such as bars, rings, etc, must be accounted for to improve our knowl¬ 
edge on bulge formation and evolution. 

The study of the intrinsic shape of elliptical galaxies at high redshift has recently 
suffered a boost t hanks to the arrival of high spatial resolution surveys on large 
fields of view (see Chang et'dl 20131 and references therein). This kind of stud¬ 
ies can provide an in-situ view of galaxy evolution and their application to the 
intrinsic shape of bulges will be key to further progress on this topic. 
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